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Abstract. Fast thermalization and a strong build up of elliptic flow of QCD matter were investigated within 
the pQCD based 3+1 dimensional parton transport model BAMPS including bremsstrahlung 2 «-> 3 processes. 
Within the same framework quenching of gluonic jets in Au+Au collisions at RHIC can be understood. The 
development of conical structure by gluonic jets is investigated in a static box for the regimes of small and large 
dissipation. Furthermore we demonstrate two different approaches to extract the shear viscosity coefficient from 
a microscopical picture. 



1 Introduction 

The large value of the elliptic flow U2 measured in experi- 
ments at the Relativisitc Heavy Ion collider (RHIOf 1 I2I3II 
suggests, that in the hot and dense fireball equilibration of 
quarks and gluons occurs on a very short time scale < 1 
fm/c. This results into the assumption, that this locally ther- 
malized matter, the quark gluon plasma (QGP), behaves 
like a nearly perfect fluid. In this case quarks and glu- 
ons should be rather strongly coupled and the QGP should 
have a very small viscosity to entropy ratio rj/ s. This value 
might in fact be close to the conjectured lower bound rj/s = 
1 1 (47r), obtained from a correspondence between confor- 
mal field theory (CFT) and string theory in an Anti-de- 
Sitter space (AdS) g]. 

The phenomenon of jet-quenching was another impor- 
tant discovery at RHIC Q. Hadrons with high transverse 
momenta are suppressed in Au +Au collisions with respect 
to a scaled p + p reference 06171 . This quenching of jets is 
commonly attributed to energy loss on the partonic level as 
the hard partons produced in initial interactions are bound 
to traverse the QGP created in the early stages of heavy-ion 
collisions (HIC). In addition, very exciting jet-associated 
particle correlations were observed [8 j, which might be the 
result of a conical emission off propagating shock waves 
in form of Mach Cones. These Mach Cones might be in- 
duced by high-energetic partons traversing the expanding 
medium |9|. 

A large class of phenomena in heavy-ion collisions can 
be investigated within the framework of the kinetic trans- 
port theory. Among others, the relativistic pQCD-based 
on-shell parton transport model BAMPS (Boltzmann Ap- 
proach to Multiparton Scatterings) [10] was developed to 
describe the early QGP stage of a heavy-ion collision (HIC). 
Using BAMPS early thermalization of gluons within t < 1 
fm/c was demonstrated in Au+Au collisions at y/s^N = 
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200 GeV employing Glauber initial conditions of mini- 
jets and the coupling constant as = 0.3 - 0.6. In addi- 
tion to the elastic pQCD gg <-> gg processes, pQCD-based 
bremsstrahlung gg <-> ggg was included. These were shown 
to be essential for the achievement of local thermal equilib- 
rium at that short time scale. The fast thermalization hap- 
pens also in a similar way with Color Glass Condensate 
initial conditions [ 1 1 1. 

BAMPS has been recently applied to calculate elliptic 
flow and jet quenching at RHIC energies |12l for the first 
time in a consistent and fully pQCD-based microscopic 
transport model. Both key observables could be addressed 
on the partonic level within a common setup. The left panel 
of Fig. [U shows that the medium simulated in the parton 
cascade BAMPS exhibits a sizable degree of elliptic flow 
in agreement with experimental findings at RHIC as dis- 
cussed in Ref. 013I14L 

Extraction of the shear viscosity over entropy density 
ratio J]/ s from BAMPS simulations confirmed the essen- 
tial importance of inelastic processes. Bremsstrahlung and 
back reaction processes lower the shear viscosity to en- 
tropy density ratio significantly, by a factor of 7, compared 
to the values obtained if only elastic collisions are consid- 
ered [JJ 16|. For Of, = 0.3 one finds rj/s - 0.13, whereas 
for as - 0.6 the obtained value matches the lower bound 
of 77/s - l/(47r) from the AdS/CFT conjecture. 

In these proceedings we show some recent descriptions 
and developments of different phenomena in relativistic 
HIC using BAMPS. Due to the large momentum scales in- 
volved, the energy loss of partonic jets can be treated in 
terms of perturbative QCD (pQCD) and most theoretical 
schemes attribute the main contribution to partonic energy 
loss to radiative processes |17|. In addition, the possible 
propagation of Mach Cones in the QGP induced by such 
high-energetic partons can be studied. Considering the ear- 
lier works, in which the effects of dissipation on relativis- 
tic shock waves were investigated 1181191201 . we demon- 
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stxate the transition of Mach Cones from ideal to the vis- 
cous ones. It is a major challenge to combine jet physics on 
the one hand and bulk evolution on the other hand within 
a common framework. In the end we demonstrate two new 
independent methods to extract numerically the shear vis- 
cosity to entropy density ratio rfl s. One of them is based 
on the classical picture of a linear velocity gradient, where 
the shear viscosity coefficient 77 can be calculated from the 
local deformation of the distribution function. The other 
method is based on the Green-Kubo relation, where shear 
viscosity is obtained from the autocorrelation function of 
equilibrium fluctuations. 



up to roughly pj ~ lOGeV, see left panel of Fig. [2] For 
this a critical energy density Sc - 0.6GeV/fm^ has been 
used, in order to be comparable to previous calculations of 
the elliptic flow within BAMPS. The V2 of high-/?^ gluons 
is at first rising with pj, but from pj ^ 4 to 5 GeV on, it 
begins to slightly decreases again. This behavior is in good 
qualitative agreement with recent RHIC data |23 | that for 
charged hadrons shows V2 to be rising up to V2 ~ 0.15 at 
Pt ~ 3 GeV followed by a slight decrease. 

3 Transition from ideai to dissipative i\/lacli 
Cones 



2 Jet Quenchiing in Au+Au coliisions at 200 
AGeV 

For simulations of Jet Quenching in heavy ion collisions 
the initial gluon distributions are sampled according to a 
mini-jet model with a lower momentum cut-off po = lA 
GeV and a /T-factor of 2. The test particle method IfTOl 
is employed to ensure sufficient statistics. Quarks are dis- 
carded after sampling the initial parton distribution since 
currently a purely gluonic medium is considered. To model 
the freeze out of the simulated fireball, free streaming is ap- 
plied to regions where the local energy density has dropped 
below a critical energy density Sc (se = l.OGeV/fm"' un- 
less otherwise noted). 

The right panel of Fig. [1] shows the gluonic Raa sim- 
ulated in BAMPS for central, b - Ofm, collisions. It is 
roughly constant at ~ 0.053 and in reasonable agree- 

ment with analytic results for the gluonic contribution to 
the nuclear modification factor Raa I 17|, though the sup- 
pression of gluon jets in BAMPS appears to be slightly 
stronger We expect improved agreement in future studies 
when employing a carefully averaged {b) that will be better 
suited for comparison to experimental data than the strict 
= fm case. 

We have computed the gluonic Raa for non-central Au 
+ Au collisions at the RHIC energy of ^fs - 200AGeV 
with a fixed impact parameter = 7 fm (Fig. |2]i, which 
roughly corresponds to (20 - 30)% experimental central- 
ity. A comparison in terms of the magnitude of the jet 
suppression for b = 7 fm is difficult since there are no 
published analytic results available to compare to. Tak- 
ing the ratio of the b = 7 fm to the = fm results as 
a rough guess indicates that the decrease in quenching is 
more pronounced in BAMPS compared to experimental 
data. The ratio of the nuclear modification factor between 
central (0-10)% and more peripheral (20-30)% collisions 
is RAA\im-io% I ^/1aI2o%-3{)% ~ 0.6 for the experimental 
data, while for the BAMPS results RAA\h=oim I RAA\b=ifm ~ 
0.4. However, the issue of detailed quantitative comparison 
needs to be re-investigated once light quarks and a frag- 
mentation scheme are included into the simulations. 

To complement the investigations of Raa at a non-zero 
impact parameter b - 1 GeV, we have computed the el- 
liptic flow parameter V2 for gluons at the same impact pa- 
rameter and extended the range in transverse momentum 



In the early 1970s shock waves were theoretically predicted 
to occur in relativistic heavy-ion collisions (HIC) ||24j and 
they have been experimentally investigated. There is an im- 
portant issue in recent studies whether the small but finite 
viscosity allows the development of relativistic shocks in 
form of Mach Cones in such a hot and dense matter like the 
QGP. Within the framework of BAMPS studies were fin- 
ished to answer the question whether shocks can develop 
with finite viscosity and how this will alter such a pic- 
ture 1 18J . Within the relativistic Riemann problem it was 
shown that one dimensional shocks smears out if viscosity 
is large ri9'201. However, the expected viscosity in HIC 
seems to be small enough to allow a significant contribu- 
tion of shocks in form of Mach Cones. In the following we 
report a very recent study. 

Mach Cones, which are special phenomena of shock 
waves, have their origin in ideal hydrodynamics. A very 
weak perturbation in a perfect fluid induces sound waves 
which propagate with the speed of sound c, = ^[dpjde, 
where p is the pressure and e is the energy density. In the 
case where the perturbation with velocity ujet propagates 
faster than the generated sound waves, the sound waves lie 
on a cone. Considering a gas of massless particles, where 
e - 3p and Cj = 1 / V3, then the emission angle of the 
Mach Cone is given by a^, = arccos(cj/t;jet) = 54, 73°. 

A stronger perturbation induces the propagation of shock 
waves exceeding the speed of sound, therefore the emis- 
sion angle changes and can be approximated by 
a ^ arccos((;shock/fjet)- Here 



^shock — 



iPi - PQ)(eo+pi) 



(ei - eo)(ei + po) 



is the velocity of the shock front, po (eo) the pressure (en- 
ergy density) in the shock front region and pi (ei) in the 
stationary medium itself. The expression has the following 
limits: For po » pi we obtain Ushock ~ 1, whereas for a 
small perturbation, po ~ pi, we get Ushock ~ Cs- 

We employ the microscopic transport model BAMPS 
to investigate Mach Cones with diff'erent strength of dis- 
sipations in the medium using a jet moving in positive z- 
direction, initialized at f = fm/c at the position z - -0.8 
fm. The jet is treated as a massless particle with zero spatial 
volume and zero transverse momentum, that is, p. - Zijet = 
200 GeV and yjet - 1 . The energy and momentum deposi- 
tion to the medium is realized via collisions with medium 
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Fig. 1. Left panel: Elliptic flow V2 as a function of the number of participants for Au+Au at 200 AGeV for different combinations of the 
strong coupling ff, and the critical energy density e^. See (T4| for more information. 

Right panel: Gluonic Rj!^ at midrapidity (j/ e [-0.5, 0.5]) as extracted from simulations for central Au+Au collisions at 200 AGeV. For 
comparison the result from Wicks et al. 0171 for the gluonic contribution to Raa and experimental results from PHENIX 11211 for n° and 
STAR |5 1 for charged hadrons are shown. 
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Fig. 2. Left panel: Elliptic flow V2 for gluons in simulated Au+Au collisions at 200 AGeV with b = 7fm. e^, - 0.6GeV/fm^. See f22\ 
fore infoiTnation. 

Right panel: Gluonic Raa as extracted from BAMPS simulations for & = Ofm and b = 7fm, = 1.0GeV/fm\ For comparison 
experimental results from PHENIX [21] for are shown for central (0 - 10)% and off-central (20 - 30)% collisions. 



particles. In this scenario we neglect the deflection of the 
jet and it can not be stopped by the medium; its energy and 
momentum is set to its initial value after every collision. 

All simulations are realized within a static and uniform 
medium of massless Boltzmann particles and T = 400 
MeV. For this study we consider only binary scattering 
processes with an isotropic cross section among the bulk 
particles. To save computational runtime we reduce our 
problem to two dimensions. Here we choose the ^z-plane 
and apply a periodic boundary condition in (/-direction. 

In Fig[3]we demonstrate the transition from ideal IVIach 
Cone to a highly viscous one by adjusting the shear viscos- 
ity over entropy density ratio in the medium from tj/s = 
l/647r 0.005 to l/n x 0.32. The energy deposition of 
the jet is approximately dE/dx = 11-14 GeV/fm. We 
show a snapshot at f = 2.5 fm/c. 

Using a small viscosity of ri/s = l/647r, we observe 
a strong collective behavior in form of a Mach Cone, as 
shown in the left panel of FigE] Due to the fact that the en- 



ergy deposition is strong, the shock propagates faster than 
the speed of sound through the medium. For comparison, 
the ideal Mach Cone caused by a very weak perturbation is 
given by a solid line with its emission angle = 54, 73°. 
Furthermore, a strong diffusion wake in direction of the 
jet, characterized by decreased energy density, and a head 
shock in the front are clearly visible. 

If we increase the viscosity of the medium to larger 
values, shown in the middle and left panel of Fig[3] the 
typical Mach Cone structure smears out and vanishes com- 
pletely. Due to stronger dissipation, the collective behav- 
ior gets weaker because of less particle interactions in the 
medium with a larger tj/s. The results agree qualitatively 
with earlier studies 11181201 . where a smearing-out of the 
shock profile is observed with higher viscosity. 

In addition to the scenario of a punch through jet we 
demonstrate in Fig|4] the scenario of a massless jet with 
finite energy which can also be deflected. Its initial energy 
is set to p, = Zijet - 20 GeV, where the starting point is 



Tj/s = 1/64JC 
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Fig. 3. (Color online) Scenario of a massless jet with p- = Ej^ = 200 GeV in the punch through scenario - the shape of a Mach Cone 
shown for different viscosities of the medium, r]/s = l/64n (left), rj/s = l/4n (middle), rj/s = l/n (right). The energy deposition of the 
jet is approximately dE/dx = 11-14 GeV/fm. We show the energy density plotted together with the velocity profile. Additionally, in 
the left panel the linear ideal Mach Cone for a very weak perturbation is shown by a solid line; its emission angle is a,„ = 54, 73°. 



Ti/s = 0.5 




Fig. 4. (Color online) Scenarion of a deflectable jet with finite energy p, = £jet = 20 GeV - Induced Mach Cone structure for different 



viscosities of the medium, 
velocity profile. 



rj/s = 0.01 (left), ri/s = 0.1 (middle), rj/s = 0.5 (right). We show the energy density plotted together with the 



z - -0.3 fm. We show the results for different viscosities, 
Tj/s X 0.01 to 0.5 att- 5.0 fm/c. In analogy to the results 
above we observe a clear Mach Cone structure for small 
viscosities and a smearing out with larger values of rj/s. 
Only in the ideal case a strongly curved structure in which 
the building up of a strong vortex is visible. The physical 
meaning of these phenomena and also jets with the full 
pQCD cascade have to be explored in future studies l25l . 

4 Extraction of shear viscosity from 
microscopic tlieory 

Whereas in viscous hydrodynamics the shear viscosity to 
entropy density ratio rj/s is the important input parameter, 



in microscopic transport theory the cross section or colli- 
sion rate governs the behavior of the system. A big chal- 
lenge is to connect these both parameters e.g. by extracting 
the shear viscosity rj from transport calculations. There is a 
big interest to calculate the shear viscosity in a full pQCD 
simulation including inelastic processes, which was done 
in previous calculations II15I16I . In this section we present 
another two different methods to extract the shear viscosity 
T] numerically from transport calculations employing the 
transport model BAMPS. 

The first method is motivated by the classical picture 
of a stationary velocity gradient. We introduce a particle 
system embedded between two plates moving in opposite 
z-direction represented by thermal reservoirs of particles 
with +y„aii. In x-direction the system has a size of L, in 
y and z direction the system is homogeneous. We assume 
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T [fm/c] time [fm/c] 

Fig. 6. Left panel: The correlation function for different isotropic cross section is shown. The relaxation time t is linear dependent on the 
cross section a. 

Right panel: A snapshot of a fluctuating shear stress tensor, which is used to calculate (in this case for pQCD interactions) the correlation 
function. 



the mean free path of the particles to be small compared to 
the system size, i.e. /l^fp << After a sufficient long time 
scale a stationary velocity gradient v^ix) is established. 

In the non-relativistic limit the velocity gradient is lin- 
ear, dv2,{x)ldx - const. Within the Navier-Stokes Ansatz 
the shear stress tensor in x-z- direction is proportional to 
the gradient of the velocity: 

dx 

In the true relativistic limit the velocity gradient is not lin- 
ear, but has the shape of a hyperbolic tangent 

v{x) - tanh(6'jic) 

where the rapidity ijzix) - 9x is the argument of the hyper- 
bolic tangent and 9 - dy^ldx. 

We employ BAMPS to extract the velocity profile with 
fwaii = 0.964. We use a massless Boltzmann gas, where 
the interactions between particles are realized only via bi- 
nary collisions and isotropic cross sections. The results are 
shown in Fig|5] On the left panel we show the velocity pro- 
file, which has the expected shape of a hyperbolic tangent. 



The results are compared to a linear fit with tanh(0.4x) ac- 
cording to the expected rapidity profile shown on the right 
panel of Fig|5] For small velocities the system achieves the 
correct non-relativistic limit, where the velocity gradient is 
linear 

Extracting the gradient of the rapidity profile 6 allows 
us to calculate the shear viscosity coefficient using the com- 
mon relativistic relation for the shear stress tensor 

n-'"- = -riye, 

where tt'-' for a kinetic cascade can be extracted numeri- 
cally (right hand side) from BAMPS 

d'p r)=y ^5(rj - r) 

P "pi P j 

To verify the accuracy of this method the results were com- 
pared to the analytical expression f26l for binary collisions 
and isotropic cross sections 77 = 1.2654 ■ T/cr, where T is 
the temperature and cr is the implemented cross section. 
We found an excellent agreement. 

The use of Green-Kubo relations is another possibility 
to extract the shear viscosity from numerical simulations. 
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collision type 


{r]l s) shear flow 


(77/1) green kubo 


(77/i)Xuetal. fTSl 


(77/ .s) El etal. |16 | 






1.035 + 0.015 


1.03 


1.4 


2 2 and 2 «-> 3 


0.12 ± 0.07 


0.13 ± 0.02 


0.13 


0.16 



Table 1. rjl s values extracted from BAMPS using the shear flow and Green-Kubo relation discussed in this proceeding. Full pQCD 
simulations with and without inelastic 2 «-» 3 Bremsstrahlung processes were performed with or, = 0.3 and T = 0.4 GeV. The extracted 
values are compared to the results of previous calculations discussed in I15I16I . 



Green-Kubo relations connect a linear transport coefficient 
with the integrated correlation function of the underlying 
flow connected with the transport coefficient. In hydrody- 
namics the shear stress tensor is generated by a shear gradi- 
ent. Using the Green-Kubo framework the shear viscosity 
7] is found to be linear with the correlation function of the 
equilibrium shear tensor The presence of fluctuations is at- 
tested to every equilibrated system and Onsager's regres- 
sion hypothesis |27| tells that fluctuations are driven back 
to equilibrium by the same transport parameter as small 
non-equilibrium deviations, which motivates the Green- 
Kubo relations. The type of interaction and its property, 
for example the collision rate, directly influence the emer- 
gence and decline of fluctuations. Therefore these effects 
also influence the shape of the correlation function which 
is flie link to the shear viscosity. With a semi-analytic dis- 
cussion of the correlation function the shear viscosity of 
the system can be reduced to be only dependent on the 
mediums relaxation time of its auto correlation function. 
As a consequence, this relaxation time includes all interac- 
tion and medium effects. 

In relativistic notation the Green-Kubo relation for the 
shear viscosity has the following form [28 1: 

1 V r°° r r 

n = lim — ■ -— — ; I dT I dr I dx 

X J dt {n'-ix, t) ■ TT^-ix + r,t + T)) , 

where 7r*~(r, t) is defined above. 

The expression for the shear viscosity can be brought 
in a more convenient form when the correlation function 
of shear tensor is known. BAMPS solves the Boltzmann 
equation for an ultra relativistic gluon gas with stochas- 
tically interpreted cross sections. The observed correlation 
function follows an exponential decay Cor(T) = Co-exp-f/r. 
The zero time correlation Co can be derived analytically 
by the variance of the equilibrium Boltzmann distribution, 
the equilibrium gluon density and its projection on the xy- 
plane: 

Cor{T) = ^ Cor{r, t = 0)dr = ^- j ■ ■ exp -f/r 

With the exponential ansatz and the zero time coiTelation, 
the Green Kubo relation can be integrated analytically. To 
calculate the shear viscosity to entropy density, the equi- 
librium density for an ultra relativistic gas is used: s = 
An. This will give the simplified Green-Kubo relation for 
BAMPS, which depends only on t: 




T includes all interaction effects of the medium, see Fig. 
|6] for a correlation functions with isotropic and constant 
cross sections. Again, this method is tested for binary and 
isotropic cross sections and comparison gives an excellent 
agreement to the known results. 

Finally, in Table [1] we now summarize the results ex- 
tracted from full pQCD simulations performed with BAMPS 
using both introduced methods. We use a coupling con- 
stant ffj - 0.3 and T - 0.4 GeV and show the results 
with and without inelastic pQCD Bremsstrahlung. In ad- 
dition the results are compared to previous calculations, 
based on Navier Stokes ||15J and Grads method [16|. The 
results presented in lfT6l are extracted from full dynami- 
cal simulations, therefore loss of chemical equilibrium and 
the size of equilibrium deviations might have an effect on 
the obtained values. In contrast, the methods discussed in 
these proceedings employ a static setup. We observe an 
overall good agreement. For a, - 0.3 the shear viscosity 
to entropy density ratio is around 77/s = 0.13 as raised in 
the introduction and in |15|; this low value is mainly due 
to the implementation of inelastic 2 «-» 3 Bremsstrahlung 
processes |15|, which are thus of crucial relevance to un- 
derstand the perfect liquid behavior of the system. 

A future challenge will be to extract the shear viscosity 
77 for much small coupling constant a, and to check ||29l 
for the expected rj ~ a] log(l/ai) behavior Ii30j . 
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